Intervalley splitting and intersubband transitions in n-type Si/SiGe quantum wells: 

pseudopotential vs. effective mass calculation 
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Intervalley mixing between conduction band states in low-dimensional Si/SiGe heterostructures 
induces splitting between nominally degenerate energy levels. The symmetric double-valley effective 
mass approximation (DVEMA) and the empirical pseudopotential method (EPM) are used to find 
the electronic states in different types of quantum wells. A reasonably good agreement between 
the two methods is found, with the former being much faster computationally. Aside from being 
an oscillatory function of well width, the splitting is found to be almost independent of in-plane 
wave vector, and an increasing function of the magnitude of interface gradient. Whilst the model 
is defined for symmetric envelope potentials, it is shown to remain reasonably accurate for slightly 
asymmetric structures such as a double quantum well, making it acceptable for simulation of mul- 
tilayer intersubband optical devices. Intersubband optical transitions are investigated under both 
approximations and it is shown that in most cases valley splitting causes linewidth broadening, 
although under extreme conditions, transition line doublets may result. 

PACS numbers: 73.21. Fg, 78.67.De 
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I. THEORETICAL CONSIDERATIONS 

A. Intervalley Mixing and the Effective Mass 
Approximation 

The conventional single- valley effective mass/envelope 
function approximation (EMA) is computationally effi- 
cient and remains remarkably accurate at the atomic 
scale. 1 The behaviour of carriers in a 2- dimensional het- 
erostructure is modeled by introducing an effective mass, 
m*(z) and an envelope function potential, V{z) which 
depends on the material composition-. The band non- 
parabolicity is optionally incorporated into the effective 
mass, via m*(z) = mo [1 + a(E — V(z))]. The Hamilto- 
nian for such a system is: 
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If there are multiple equivalent valleys in fc-space, these 
must be considered independently within the "pure" 
EMA, with any intervalley mixing effects neglected. The 
conduction band edge is located in the six A-valleys for 
Sii-zGe-E alloys with x < 85%2. Strain effects split 
the degenerate valleys into four Ay -valleys and two A^- 
valleys. As the A-valleys are ellipsoidal, a different effec- 
tive mass is used for each degenerate set. 

With two equivalent A_L-valleys, their quantum con- 
finement subbands are predicted by the envelope func- 
tion/EMA to be degenerate. In reality, there is a coupling 
between the two sets of states, which results in a split- 
ting, i.e. lifting of the degeneracy. This effect has been 
observed experimentally in Shubnikov-de Haas oscillation 
measurements in high magnetic fields^&Z ' 8 ' 9 ' 10 ' 11 with 
energy splitting up to a few meV. Boykin et al presented 



a tight-binding model of the ground state splitting in a bi- 
ased square quantum well with both hard-wall and cyclic 
boundary condition o 12 : 13 ' 14 . The ground state splitting 
in an unbiased square well was shown to be approximated 
by: 



16tt 2 u 
(5 + 2)3 



■ sin 



sm 



(5 + 2)- 



(2) 



where (j)mm = kga, and fco denotes the position of the 
valley minimum in the Brillouin zone, a is the lattice 
constant, S is the number of crystalline monolayers in 
the quantum well and u is a fitting constant. From this 
model, it is clear that the ground state splitting oscil- 
lates with well width; the frequency being dependent on 
the location of the valley minima. Similar results have 
been obtained for the two lowest subbands in an unbi- 
ased well by Chiang^ in an anti-bonding orbital model, 
and recently by Nestokloni^ in a slightly different tight- 
binding model. 

Splitting due to an electric field has been considered 
by modelling a triangular quantum well. The compo- 
sition profile of the structure in this case is selected to 
provide a potential gradient i.e. an internal electric field. 
Although an effective mass model by Shami 7 . proposed 
that the splitting is simply proportional to the applied 
field, Boykin et ai^ and Grosso et al— show that the 
splitting is a non-linear function of both the well width 
and the electric field. 

Modifications to the EMA have been introduced in or- 
der to describe the intervalley mixing effects in an infinite 
square wel l 20 ' 21 : 22 and a finite square well with impu- 
rity states^ 3 -. Ting and Chang's Double Valley Effective 
Mass Approximation (DVEMA) 2 - provides an elegant 
self-contained EMA description of intervalley mixing ef- 
fects. This model allows a computationally fast Simula- 
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tion of long structures such as Quantum Cascade Lasers 
(QCLs)2£, where the use of atomistic methods would 
be extremely cumbersome. The model is applicable to 
any symmetric conduction band-edge envelope potential, 
V(z). The kinetic energy operator in the Hamiltonian is 
unchanged, but a splitting potential U (z) is included to 
account for intervalley mixing: 



h 2 d 



1 d 



2 dz \m(z) dz 



+ V(z)±U(z). (3) 



The derivation of Ting and Chang's DVEMA^i starts 
with the assertion that in fc-space there are two equiva- 
lent A_l -valleys at wave- vectors, k ± k$. By symmetry, 
the wavefunction can be decomposed into either even or 
odd symmetric combinations of functions centred at the 
valley minima with equal and real weighting coefficients: 
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In the most general case, the weighting factors of the two 
terms may become complex conjugates, with magnitude 
) 26 . The complete wavefunction is now defined as: 

v2 



|V>i,a> = J2 Mk)\h, 



(5) 



The matrix elements of the envelope potential operator 
in the Hamiltonian are written in the basis defined by 
Eq. Q: 

V nm = (k n \V \k m ) = \ ((k n + k \±(k n - k \) 

2 (6) 

x V (\k m + k ) ± \k m - k Q )) 
Rearranging this expression yields: 

V nm = ~ ((k n + k Q | V \k m + fc ) + (k n ~ k \ V \k m - fe )) 

± - ((k n + k \ V \k m - ho) + (k n -k \V \k m + k )) 

(7) 

Using the discrete to continuous approximation, V nm — 
V(k m — k n ), the intervalley envelope term can be written 



Vi, 2 (k) = V(k)±-\V(k- 2k Q ) + V(k + 2k ) 



(8) 



where V(k) is the Fourier transform of the conduction 
band-edge envelope potential. The configuration-space 
form of the intervalley envelope function is found by tak- 
ing the inverse Fourier transform of this result and the 
splitting potential is therefore extracted as: 



U(z) 
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where T~ x denotes the inverse Fourier transform. This 
form of real-space splitting potential is more general than 
the special case of a square well with abrupt interfaces^. 
It applies to any symmetric potential, but does not ex- 
plicitly show that the splitting potential oscillates as a 
function of well width W of the form sin(fcoWO/fc(j 24 . 
In addition, the numerical solution is complicated by 
the presence of a continuous, rapidly varying splitting 
function, as opposed to a pair of delta-functions at the 
square well interfaces 2 ^. However, since the problem is 
one-dimensional, this is not a major concern. For asym- 
metric envelope functions, the simple cosine form is no 
longer strictly valid as the Fourier coefficients in equation 
2] become a complex conjugate pair. However, it will be 
shown here that the symmetric approximation still gives 
excellent results for structures with a moderate degree of 
asymmetry. 



B. Empirical Pseudopotential Calculation 

An empirical pseudopotential method (EPM) is used 
here to calculate electronic states in Si/SiGe based quan- 
tum well structures, and provide a comparison with the 
results of the DVEMA simulation. Being a microscopic 
method (in common with tight- binding) , it normally re- 
veals intervalley interference induced splitting of size- 
quantised subbands, without any adjustable parameters 
introduced for this purpose. The supercell implemen- 
tation of the EPM was used, with a continuous atomic 
formfunction, V(g). The "modified Falicov" formfunc- 
tion described by Friedel et al^ was selected: 
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This formfunction has been used by Fischetti and Laux2&, 
as well as in previous work by the authors 2 ^. It gives rea- 
sonable agreement with experimental data for both bulk 
Si and Ge band structure and for band discontinuities at 
the interface. A cut-off energy of 4.5 Ry was used, which 
gives an acceptable number of plane waves for accurate 
and rapid computation with all the structures considered. 
The parameters for Si and Ge are given in Table HI and 
the virtual crystal approximation was used for the alloy. 

The EPM can be used for structures with either abrupt 
interfaces or graded compositions. In the latter case, 
the interface grading is piecewise constant (i.e. within 
the crystalline monolayer width). Individual layers are 
given the required alloy compositions. In contrast, the 
DVEMA uses a continuous potential profile, though fea- 
tures smaller than the width of a monolayer have little 
practical meaning. 

It is important to note that effective mass based cal- 
culations (like the DVEMA) can never fully reproduce 
the results of microscopic EPM modelling. This is be- 
cause the DVEMA only handles four bulk states explic- 
itly, while the EPM implicitly includes many evanescent 
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TABLE I: Pseudopotential parameters . 



Parameter 


ai 


ai 


03 


ci4 


0-5 


(Z6 


Si 


212.1372 


2.2278 


0.6060 


-1.9720 


5.0 


0.3 


Ge 


108.9024 


2.3592 


0.7400 


-0.3800 


5.0 


0.3 



states, stemming from remote bulk bands. Furthermore, 
the location of indirect valleys may vary between the bulk 
materials used in the well and barrier — a situation which 
is difficult to handle with effective mass methods. Given 
that remote bands are usually less important than the 
bands from which the quantised states are derived, one 
can expect reasonable accuracy from the DVEMA. This 
may be validated by comparison against the EPM calcu- 
lation. 



TABLE II: Material parameters for silicon and germanium 

Parameter Si Ge Unit 

o 543.1 a 563.3° pm 

Cn 167.5 6 131.5 C GPa 

C12 65.0 6 49.4 C GPa 

"Reference 3^. 
b Reference 33. 
c Reference 33- 



II. NUMERICAL RESULTS AND DISCUSSION 

DVEMA and EPM calculations were performed for a 
range of Si/SiGe Quantum Wells (QW), using the ma- 
terial parameters as follows. The lattice constant a{x) 
of relaxed Si^Gei-^, alloy is found by interpolation of 
the elemental lattice constants, given in table |TT] accord- 
ing to a(x) — asi(l — x) + acoX — fobow^(l — x)^, with 
the bowing parameter &bow = 0.2733 pm2£. In a struc- 
ture coherently grown in the [001] direction, compris- 
ing layers of different compositions on a relaxed virtual 
substrate, the strain results in the perpendicular lattice 
constant of a particular layer being equal to a±(x) — 
a(x) [1 - (2C 12 (x)/C n (x)) ((a (x) ~ a(x))/a{x))], while 
the in-plane lattice constants equal that of the substrate. 
The elastic constants for SiGe are found by linear inter- 
polation of the values given in Table ITTl and ag(x) is the 
bulk lattice constant of the relaxed virtual substrate on 
which the structure is grown. 

In this work, the Ge fraction in the virtual substrate is 
fixed at 20%. For the DVEMA calculation, the A- valley 
longitudinal and transverse effective masses are uil = 
0.916 and vht = 0.19 respectively, and do not depend 
on the alloy composition. In the EPM calculations the 
total length of the structure (i.e. the supercell period) 
which includes the well and barrier layers, was set to a 
fixed value of 35 monolayers (ML) , where 1 ML is half 
the lattice constant (i.e. twice the crystalline monolayer 
width of eqn. [5]). This maintained a constant number of 
plane- waves in the pseudopotential basis set, and avoided 
fluctuations in the results caused by variable size of basis. 
This is important since intervalley splitting is relatively 
small on the energy scale covered by EPM^I. 



A. Finite Square Well 

The first set of calculations was for a simple square 
QW with abrupt interfaces. Fig. [1] shows the influence 
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FIG. 1: Intervalley splitting, AE in the lowest two subbands 
of an 8 ML quantum well as a function of barrier composition. 



of barrier composition (potential height) on the split- 
ting of the lowest two subbands, obtained by both the 
DVEMA and EPM calculations, for a fixed, 8 ML wide 
Si quantum well. The confining potential of the quan- 
tum well increases almost linearly with the Ge content in 
the barriers^. The results show a clear increase in valley 
splitting with increasing confining potential, with the two 
sets of results being broadly in good agreement, and with 
the most obvious discrepancy being the discontinuities in 
the DVEMA plot. 

The effect of well width upon state splitting was in- 
vestigated next, for a structure with a fixed barrier com- 
position of 50% Ge. The well width was varied between 
1 ML and 25 ML in a supercell of total length 35 ML. 
This leaves a minimum 10 ML barrier region, thus en- 
suring decoupling of neighbouring quantum wells under 
periodic boundary conditions. It also represents the re- 
alistic range of well widths within intersubband optical 
devices. Figure [H shows the EPM and DVEMA results. 
As predicted by equation^ the valley splitting is a de- 
caying oscillatory function of well width, originating from 
interference of the wavefunction components reflecting at 
the quantum well interfaces. Since the number of mono- 
layers in the structure is restricted to the set of positive 



4 





(EPM) 




(EPM) 


«-* AK 


(DVEMA) 


*-*AE 2 


(DVEMA) 




5 10 15 

Well Width (ML) 



20 



25 






(EPM) 




(EPM) 


«-* AF, 


(DVEMA) 


^* AE 2 


(DVEMA) 



10 15 
Well Width (ML) 



20 



FIG. 2: Intervalley splitting in a Si/Sio.sGeo.s quantum well 
as a function of well width 

integers, the splitting function is undersampled to show 
the precise value of the period of oscillations. DVEMA 
shows good agreement with the EPM results for the enve- 
lope of the splitting, whilst the oscillatory component is 
approximately correct. For very small well widths how- 
ever, the results for higher subbands deviate from the 
theory, as the states are no longer bound in the well. 

The tight-binding results presented by Boykin 12 con- 
sider wide quantum wells, with low-energy ground state 
and effectively infinite barriers. The decaying oscilla- 
tory form is however, the same as that obtained from 
DVEMA calculations. Equation [2j extracted from the 
tight-binding model, predicts a period of around 6 ML as 
observed in fig. [2] Setting u ss 3 yields a match in the 
amplitude between the three models. This figure is some- 
what higher than the value given in the reference above, 
although equation [2] has been fitted to the DVEMA re- 
sults for loosely bound states in a finite square well as 
opposed to being derived from bulk dispersion character- 
istics. 



B. The influence of in-plane Wave vector 

For nonzero in-plane wave vector fen , the size-quantised 
subbands in QWs generally acquire their fcy -dependent 
contributions from remote bulk bands. For bulk ma- 
terials this is manifested as non-parabolicity of the in- 
plane dispersion, and in the problem under consideration, 
fen -dependent splitting would result. The EPM calcula- 
tion automatically accounts for these effects, while the 
DVEMA would require kn -dependent "correction" terms 
as remote bands are not explicitly included in the model. 
The in-plane non-parabolicity of the A j_ -valley is rela- 
tively low, and indeed a rather weak dependence of split- 
ting on fen is found. In the EPM calculation for a 10 ML 
quantum well, the splitting in the first and second sub- 
bands increased approximately linearly by 17% and 8% 
respectively, when h\ changed from zero to 10% of the 
Brillouin zone edge (i.e. in the range with non-negligible 
electron occupancy at any reasonable temperature). This 
implies that fcy -dependent corrections in the DVEMA are 



FIG. 3: Valley splitting in lowest two subbands as a function 
of well width, W in a QW with 3-step linear graded interfaces. 
The inset shows the general structure of a linear graded QW. 

not mandatory. 

C. Graded Barrier Potential 

In real Si/SiGe quantum wells, surface segregation ef- 
fects are well documented 35 . This refers to the "prefer- 
ence" of Ge atoms to exist on the surface of the material 
rather than in the bulk during molecular beam epitaxial 
(MBE) growth, leading to a decrease in the magnitude 
of the Ge composition gradient at the nominal interfaces. 
Recent work 3 -^ has shown that atomic hydrogen etching 
reduces the effects of surface segregation, but this is im- 
practical in a multilayer structure as it is extremely time 
consuming and surface defects caused by the etching are 
likely to accumulate. It is therefore unrealistic to model 
a Si/SiGe quantum well as having abrupt interfaces. The 
effect of graded interfaces on subband splitting is there- 
fore considered. The linear-graded structure shown in 
the inset of Fig. [3] is modelled first. 

Within the EPM calculation, the linear-graded in- 
terfaces on either side of the QW are modelled as 3- 
step piecewise-linear, i.e. the interfaces spread across 
three ML, with Ge content of 17%, 33% and 50% sequen- 
tially. The results are shown in figure [3] The well width 
is defined as the full-width at half-maximum (FWHM) of 
the envelope potential. The results of DVEMA and EPM 
are in good agreement for larger well widths — those 
which allow for more than a single bound state. The 
plots show that the oscillatory component of the valley 
splitting is unchanged, although the envelope decreases 
in magnitude. This is because graded interfaces have re- 
duced large- wave- vector Fourier components in the enve- 
lope potential, which mix the two Aj_ valleys and hence 
the splitting is generally smaller. 

As the width of the graded interfaces increases, the 
splitting is further reduced, as shown in Fig.0]for a 4-step 
graded interface with Ge content of 13%, 25%, 38% and 
50%, sequentially. Again, there is a very good agreement 
between the two models. 

A linear graded interface is a somewhat idealised model 
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FIG. 4: Valley splitting in lowest two subbands as a function 
of well width in a quantum well with 4-step linear graded 
interfaces. 
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FIG. 5: Valley splitting of the lowest two subbands as a func- 
tion of well width in a quantum well with 3-step non-linear 
graded interfaces. 



as experimental evidence shows that the interface profile 
is decidedly non-linear. A 3-step grading with germa- 
nium concentrations of 13%, 38% and 50% is therefore 
used as an approximation to a typical interface compo- 
sition. The results of EPM and DVEMA calculation are 
shown in Fig. O 

The magnitude of the splitting is somewhat larger than 
for the case of linear grading, apparently because the po- 
tential gradient at the interface is now larger over a wide 
range of energies, thus corresponding to a steeper linear- 
grading at the energies of the first and second subband 
minima. The DVEMA results are again in close agree- 
ment with the EPM results. 



D. Double Quantum Well 



AE, (EPM) 



AE, (DVEMA), z =0 
*AE, (DVEMA), z =7i/(4k () ) 




5 10 15 20 

Well Width (ML) 

FIG. 6: Valley splitting in lowest subband as a function of sec- 
ond well width in the double quantum well structure, shown 
in the inset. The results are shown for EPM and DVEMA. In 
the case of the DVEMA, results are shown for two different 
origins for the symmetric approximation. 



is no longer strictly applicable. The simulated structure, 
shown in the inset of Fig. [6l has a fixed 1 ML well sep- 
arated from the second well by a 1ML, 50% Ge barrier. 
All other parameters are unchanged. 

The results for the EPM and DVEMA calculations are 
shown in Fig. [6l In this case, the structure is assumed to 
deviate only slightly from the square well, and therefore 
the periodic structure may be considered approximately 
symmetric about the z = position (i.e. the left-hand 
side of the structure shown in the inset of Fig. [6]). The 
axis of symmetry, z s is therefore set at this point. As the 
structure only contains a relatively small perturbation 
from a symmetric quantum well, the DVEMA and EPM 
results are still in good agreement. The splitting energy 
is again lower than the simple square well case, since the 
left hand side (with a thin well and a thin barrier) can 
be viewed as a "soft" , nonabrupt interface. Fig. [5] also 
shows the DVEMA results when the axis of symmetry is 
shifted to z s = such that the splitting potential be- 
comes sine-modulated as opposed to cosine-modulated. 
This represents worst-case selection of the axis of sym- 
metry, if z s = is assumed to be the best. The oscilla- 
tory component of the valley splitting now appears out 
of phase with the EPM results, although the envelope of 
the oscillations is approximately correct. The symmet- 
ric approximation is therefore dependent on the origin of 
the coordinate system. However, a good estimate of the 
range of the valley splitting is possible, even with a poor 
choice of origin. 



Next, the valley splitting in a double quantum well 
structure is considered. This relatively simple structure 
usually provides sufficient design freedom for the required 
subband spacing in an optically pumped intersubband 
laser. Such a structure is also a good test of the validity 
of the DVEMA described above, as it is asymmetric and 
the simple cosine modulated splitting envelope potential 



E. Intersubband optical transitions 

Optical matrix elements were calculated for intersub- 
band transitions in the square well (section- Ill A[) . The re- 
sults from the EPM and DVEMA simulations are shown 
in figure [7] along with the separation of the transition 
energies. The difference between the optical matrix el- 
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Well Width (ML) 



FIG. 7: Optical matrix elements as a function of well width 
for the permitted transitions between the first and second sub- 
bands in the finite square quantum well considered in section 
III Al The difference between the two permitted transition 
energies is also displayed. 



FIG. 8: Permitted optical transitions in a square quantum 
well are usually between the two upper or lower valley-split 
states (left). When close to splitting minima however, this 
situation is sometimes reversed in EPM simulation (right). 

ements is small and approaches zero as the transition 
energies converge. This implies a similar magnitude of 
spectral contribution from each pair of valley-split states. 
The two methods are in close agreement for lower well 
widths, with the DVEMA predicting larger matrix ele- 
ments at higher widths. In most cases, when considering 
valley splitting of states, the permitted optical transi- 
tions are from the upper "excited state" to the upper 
"ground state" and between the two lower states. How- 
ever, when close to the splitting minima (at well widths 
of 17 ML and 23 ML), the converse situation sometimes 
applies with the EPM (fig.EJ). The DVEMA always finds 
transitions to be of upper— > upper and lower— dower char- 
acter. Transitions exhibit linewidth broadening by in- 
terface roughness and carrier scattering, typically of the 
order 5-10 meV. In the majority of cases, valley split- 
ting is relatively small and will only cause an increase 
in linewidth broadening by the amount shown in figure 
[7J However, when the valley splitting is large (for exam- 



ple at around 8 ML well width), a transition line doublet 
may become apparent. 



III. CONCLUSION 

The DVEMA method presented by Ting and Chang2i 
has been extended to model intervalley-mixing in any 
symmetric structure. DVEMA and EPM methods have 
been used to calculate A-valley subband splitting in a 
range of symmetric and asymmetric Si/SiGe heterostruc- 
tures, with both abrupt and graded interfaces. This 
provides a much closer approximation to intersubband 
optical device structures than has been achieved previ- 
ously. The results of the two methods are in good or 
reasonably good agreement with each other as well as 
with published tight-binding results, with DVEMA de- 
manding less than 0.5% the computational run-time of 
EPM. Subband splitting of the order of 10 meV was pre- 
dicted for abrupt-interface square quantum wells in the 
range of well widths of interest for silicon intersubband 
devices£L2&21^i^2. It has been shown that this will 
typically lead to linewidth broadening, although at nar- 
row well- widths a transition line doublet may form, with 
both valley-split states contributing equally to the spec- 
trum. The effect of surface segregation was modelled by 
considering both linear and non-linear composition grad- 
ing at the interfaces. This was found to reduce the valley 
splitting, as it is dependent upon the potential gradient 
at the interfaces. Modelling of valley splitting is therefore 
important for the development of the as yet unrealised n- 
Si/SiGe quantum cascade laser. As these structures con- 
sist of many heterolayers, the extended DVEMA method 
presented in this work, is a valuable tool for rapid simu- 
lation. 

An important question for future designs of n-Si/SiGe 
intersubband lasers is whether a fast depopulation of the 
lower laser state may be assisted by scattering between 
valley-split subbands. Alternatively, by careful selection 
of well widths it may possible to design structures for 
which the valley splitting is minimized. 
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